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The thermal conductivity is calculated with the Helfand-moment method in the Lennard-Jones 
fluid near the triple point. The Helfand moment of thermal conductivity is here derived for molec- 
ular dynamics with periodic boundary conditions. Thermal conductivity is given by a generalized 
Einstein relation with this Helfand moment. We compute thermal conductivity by this new method 
and compare it with our own values obtained by the standard Green-Kubo method. The agreement 
is excellent. 
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I. INTRODUCTION 



This paper completes the analysis carried out in the 
companion paper [l| by the study of thermal conduc- 
tivity with the Helfand-moment method [2j. Since the 
seventies, several works have been carried out for the cal- 
culation of the transport coefficients. Shear viscosity is 
the most studied transport property, while thermal con- 
ductivity has been less studied. The first computation of 
thermal conductivity goes back to the work of Alder et 
al. for hard-sphere systems [3] . Soft-sphere potential sys- 
tems were considered for the first time a few years later by 
Levesque et al. [1] using the well-known Green-Kubo for- 
mula [! [3, 0, HI • Nevertheless, since the autocorrelation 
function of the microscopic flux decreases slowly, espe- 
cially near the triple point, most of the studies devoted 
to the computation of thermal conductivity were per- 
formed by nonequilibrium methods. This is usually done 
by fixing the temperature gradient and measuring the 
heat flux [1 El El El El E3 , except in Ref. [l| where 
the opposite is done, i.e., the heat flux is fixed and the 
temperature gradient is measured. However, some equi- 
librium molecular dynamics studies using the standard 
Green-kubo formula were carried out fl.Jl2.JT4l E3] ■ More 



recently, as suggested in the nineties [ljj[18|], the general- 
ized Einstein relation expressing thermal conductivity as 
the variance of the time integral of the microscopic flux 
was computed [19j . The purpose of the present paper is 
to extend the Helfand-moment method of the compan- 
ion paper [l| from viscosity to thermal conductivity. As 
an application of our method, we calculate the thermal 
conductivity of the Lennard-Jones fluid at a phase point 
near the triple point. We show that the values obtained 
by the Helfand-moment method and the Green-Kubo val- 
ues are in excellent agreement. Such a method brings not 
only some interest as an alternative equilibrium molecu- 
lar dynamics method, but also in the context of recent 
theories in nonequilibrium statistical mechanics. Indeed, 
the escape-rate formalism propose clear and direct rela- 
tionships between the transport coefficients and quanti- 
ties characterizing the microscopic chaos such as the Lya- 
punov exponents and fractal dimensions [13, HH, H3, H3] • 



Successful results have already been obtained for sim- 
ple systems: diffusion in the Lorentz gas [24j and vis- 
cosity in a two-particle system [25]. Furthermore, the 
hydrodynamic-mode method, aiming to construct at the 
microscopic level the hydrodynamic modes which are the 
solutions of the diffusion or Navier-Stokes equations also 
depend on the possibility to define Helfand moments in 
periodic systems. This construction has already been car- 
ried out for diffusion in the Lorentz gas [13] ■ By defining 
the Helfand moment, we hopefully project to extend this 
construction to the other transport processes, in partic- 
ular thermal conductivity. The paper is organized as 
follows. In Section [H] the theoretical results of Helfand 
[3] are briefly outlined. In Section IIII1 we present our 
Helfand-moment method for the calculation of thermal 
conductivity. Section [IV] gives the results of the molecu- 
lar dynamics simulations. We compare our results by the 
Heldand-moment and Green-Kubo techniques and also to 
the litterature. Finally, conclusions are drawn in Section 



II. THEORETICAL BACKGROUND 

In the fifties, the Boltzmann and Enskog theories of 
transport processes were completed by an alternative ap- 
proach developed by Green [5, 6], Kubo Q and Mori [|[. 
This new approach relates the transport coefficients to 
the time autocorrelation functions of their corresponding 
microscopic flux. In particular, the thermal conductivity 
coefficient k can be expressed as 



K = lim = — 

N,V^oc k B T 2 V 



j(«)(0)jW(i)\ dt. (1) 



where is the microscopic flux associated with ther- 
mal conductivity. This theory plays an important role 
in numerical simulations. On the other hand, Helfand 
aimed to express the transport coefficients in terms of 
generalized Einstein relations [2j and showed that ther- 
mal conductivity is given by 



lim 

Ny.t^oo 2k B T 2 Vt 



G (K) (t) -G (k) (0) 



(2) 
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where G^ (t) is the corresponding Helfand moment de- 
fined as the centroid 



N 



G^(t)=J2^a (Ea-(E a )) , 



(3) 
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of the energies of the particles 
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It can be proved that the Green-Kubo formula fl} and 
the Helfand equation ([2]) are equivalent if the microscopic 
flux is related to the Helfand moment by 



J (K) (i) 



dGW(t) 



(•5) 



Thanks to Eq. ([2]), the thermal conductivity coefficient 
is manifestly a non-negative quantity, as required by the 
second law of thermodynamics. 



Following a way similar to the case of viscosity, we con- 
clude that the Helfand moment which should be used in 
periodic systems is given by 

G w (*) = 5>. 

a 

- £I>4*> (£.-<£.}) «(i -*,) 
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The derivation is given in Appendix A. As for shear vis- 
cosity, the first extra term is due to the jumps of particles 
from a boundary to the opposite one. The second extra 
term takes into account of interactions of particles with 
the image of other particles. By adding the two extra 
terms to the original Helfand moment, we recover the 
relation G = J dtJ. 



IV. NUMERICAL RESULTS 



III. HELFAND-MOMENT METHOD 

Since the transport coefficients are bulk properties, 
they can be calculated with equilibrium molecular dy- 
namics with periodic boundary conditions. This dynam- 
ics is ruled by Newton's equations modified in order to 
take into account the periodicity: 

at m * — ' 

s 

% = £ FM ; F(r) . - , (6) 

where Arj*' is the jump of the particle a at time t s with 
|| Ari^ | = L equal to the length L of the simulation box. 
The force F(r) is here assumed of finite range r < L/2. 
The relative position appearing in the force is defined as 

r Q & = r a - r & - L a | b (7) 

where the time-dependent cell translation vector L a | h [l8| 
is chosen for the minimum-image convention \\v a b\\ < L/2 
to be satisfied. For more details, see Section III of 
the companion paper [lj. We propose here a Helfand- 
moment method for thermal conductivity in systems with 
periodic boundary conditions. As for shear viscosity 
Helfand's original expression ([3]) must be modified in or- 
der to take into account the periodicity of the system by 
the addition of extra terms: 

N 

G^(t)=^x a {E a -{E a ))+I(t). (8) 

o=l 



We carried out molecular dynamics simulation to ap- 
ply our Helfand-moment method to the calculation of 
thermal conductivity in a fluid with the standard 6-12 
Lennard- Jones potential u(r) — 4e [(cr/r) 12 — (tr/r) 6 ] . 
All calculations we performed are done with the cut- 
off r c — 2.5(7. In the reduced units with the time 
t* = ty/ e/ '(ma 2 ) and the space r* = r/a, the reduced 
thermal conductivity is given by 

a2 [™ nn\ 

K = K - — \ — . (10) 

Kb V e 

The details of the simulations are given in the preceding 
paper [l[ . To show the validity of the Helfand-moment 
method, we first compare with the results of the Green- 
Kubo formula. We depict in Fig. Q]the time derivative 
of the mean-square displacement of the Helfand moment 
(121) and the time integral of the autocorrelation function 
of the microscopic flux. The calculation is carried out 
for N = 1372 atoms at the phase point near of the triple 
point with a reduced temperature of T* = fcrjT/e = 0.722 
and a density n* = na 3 — 0.8442. As seen in Fig. [TJ the 
two methods are in perfect agreement. We estimated 
thermal conductivity by a linear fit on the mean-square 
displacement of the Helfand moment. The fit is done in 
the region between 2 and 8 time units to guarantee that 
the linear regime is reached. We calculated the thermal 
conductivity for increasing system sizes from N = 108 
atoms to N = 1372. We give the results in Table Q] and 
depict them in Fig. [2] as a function of the inverse iV _1 
of the system size. The linear extrapolation gives the 
following estimate of thermal conductivity for an infinite 
system, 

k* = 6.990 i 0.030 . (11) 



3 




o Green-Kubo 
Helfand 



0.5 



1.5 2 
time 



2.5 



3.5 



FIG. 1: Thermal conductivity at the phase point T* = 0.722 
and n* = 0.8442 for N = 1372. The plain line is the derivative 
of the mean-square displacement of the Helfand moment and 
the circles the integral of the microscopic flux autocorrelation 
function. 



We see in Table U that the values of the present work 
compare very well with the values of previous studies. 



is based on generalized Einstein relations with the 
Helfand moment defined to take into account periodic 
boundary conditions in the molecular dynamics. This 
allows us to calculate the thermal conductivity coeffi- 
cient in terms of the variance of the Helfand moment 
© ■ We show that the original Helfand moment must 
be modified by adding two extra terms. After the addi- 
tion of these extra terms, we recover the relation J = G 
between the microscopic flux J and the corresponding 
Helfand moment G. This is the generalization of the 
case of diffusion for which the velocity v x (the micro- 
scopic flux) is related to the position x (the Helfand mo- 
ment) in the same way v x — x. Consequently, we have 
showed that it is possible to calculate the thermal con- 
ductivity coefficient by the mean-square displacement of 
the Helfand moment. Indeed, our molecular dynamics 
simulations show perfect agreement between the Helfand- 
moment and Green-Kubo methods. We think that the 
new expressions for the Helfand moment given in the 
present and the companion papers [l[ solve the problems 
on the use of generalized Einstein relations in periodic 
systems. In addition to the viscosities [l[ and thermal 
conductivity, it is possible to use a similar method for 
the electric conductivity a in periodic systems. Indeed, 
by a similar derivation to the one found in Appendix A, 
one can show that the electric conductivity can be writ- 
ten as: 
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FIG. 2: Thermal conductivity at the phase point T* = 0.722 
and n* — 0.8442 as a function of the inverse of the num- 
ber N of atoms. The circles are the results of the numerical 
simulations and the dashed line the linear extrapolation. 



a = lim — — — — 

jv,y,t^oo 2k B TVt 



G {a) (t) - (G (a Ht)) 



(12) 



where the modified Helfand moment for periodic systems 
is expressed as: 



G^(t) = J2^ Z a 



Ax^6(t-t a 



(13) 



with the electric charge eZ a of the particles and 
G^(0) = 0. Finally, we notice that the Helfand-moment 
method plays an important role in the escape-rate for- 
malism |20L l2~i| and the hydrodynamic-mode method 
(2(| developed in nonequibrium statistical mechanics [12] ■ 
They establish relationships between microscopic and 
macroscopic levels in the context of the understanding of 
the origin of the irreversibility of chemical-physical phe- 
nomena. The derivation of the Helfand-moment methods 
give one the possibility to confirm numerically the theo- 
retical predictions of these formalisms. 
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APPENDIX A: DERIVATION OF THE HELFAND 
MOMENT FOR THE THERMAL 
CONDUCTIVITY IN PERIODIC SYSTEMS 

The time derivative of the modified Helfand moment 
dHJ) is given by 



dG^\t) 



dx 



dt 



dE a dl(t) 
~dT + dt 



(Al) 



where the time derivative of E a is 

dEa Pa 1 ^ 9u ab p a - p fc 



m dt 

Pa + Pb 



2 ^ dr ab 

b&a) 



m 



b(=£a) 



rn 



F(r ab ) 



(A2) 



We notice that there is here no jump in position to con- 
sider because ^f- concerns a relative position r ab = 
r a — r b — L b | a which satisfies the minimum image con- 
vention within the range of the force. Symmetrizing the 
second term in Eq. (jAlj) . we find 

dE* 1 ^ Pa+Pb 

a a^b 

1ST Pa+Pb s 
A L — ^ rr) 



a^b 

\T. L 



Pa + Pb „, , 

b F(r ab ) 



a^b 



m 

Pa + Pb 

m 



F(r ab ) (A3) 



where F(r ao ) = — F(r ba ) according to Newton's third 
law and because x a b — x a — x b — L b \ ax . Substituting this 
expression in Eq. (|A1[) , where is given by modified 
Newton's equations of motion ©, we get the expression 



r/f 



Pa+Pb , . 

6 F(r ah ) 



£$>a£) (£„-(B a »*(t-t a ) 



a^b 

l^r Pa + Pb 

a^b 



F(r a6 ) 



dt 



(A4) 



On the other hand, it is well known that the microscopic 
flux j( K ' for thermal conductivity is given by 



J {K \t) = T^(E a -(E a )) 

rn 



i a 

a^b 



b ■ -T (Tab J 



(A5) 



Since dG( ,^^ = j( K \t) by definition, we thus obtain 



that: 



dl(t) 
dt 



-^Y,^ (Ea~(E a ))5(t-t s 



J^2 L b\ax 



a^b 



Pa + Pb 



rn 



F(r ah ) 



(A6) 



Finally, the quantity to be added to the usual Helfand 
moment ([3]) is: 



a^b 



^|«— -F^). (A7) 
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Authors 


year 


method 


r cut 


T* 


N 


K* 


Ak* 


Levesque et al. \4] 


1973 


GK 


N. C. 


0.715 


256 


7.067 


0.416 










0.722 


864 


7.136 


0.277 


Evans [10] 


1982 


EM 


2.5 


0.722 


108 


6.54 


N. C. 


Massobrio and Ciccotti [12] 1984 


GK 


3.5 


0.721 


256 


6.880 


0.485 






DM 


3.5 


0.721 


256 


6.873 


0.229 


Paolini et al. [27] 


1986 


NEM 


2.5 


0.718-0.721 


average ; 


6.776 


0.166 


Heyes [14] 


1988 


GK 


N. C. 


0.72 


108 


6.7 


0.34 (5%) 












256 


6.9 


0.35 (5%) 












500 


6.5 


0.33 (5%) 


This work 


2007 HM (MD) 


2.5 


0.722 


108 


6.663 


0.10 












256 


6.867 


0.12 












500 


6.932 


0.10 












864 


6.948 


0.07 












1372 


6.946 


0.12 












oc 


6.990 


0.061 



TABLE I: Results found in the literature for the thermal con- 
ductivity in the Lennard- Jones fluid near the triple point. 
The reduced density is n* = 0.8442 except for Heyes [lit ] 
(n* = 0.848). Abbreviations: DM, nonequilibrium results 
obtained by the differential method. EM, Evans method 
flOj | . GK, Green-Kubo formula with equilibrium molecu- 
lar dynamics. NEM, results by nonequilibrium method ob- 
tained by Paolini, Ciccotti and Massobrio [27| . HM, the 
present Helfand-moment method with molecular dynamics 
(MD). N. C. means the value has not been communicated. 
a -* Thermal conductivity obtained by averaging the values for 
N = 108, 256, 500, 864. 



